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Using dimensional regularization for both infrared and ultraviolet divergences, 
we confirm that the QCD corrections to the decay width T(t — ► H + b) are equal to 
those to r(t — ► W + b) in the limit of a large t quark mass. 



Many extensions of the Standard Model contain more than the one necessary Higgs 
doublet, and the new degrees of freedom appear as extra Higgs scalars, some of which are 
charged (see ]I| and references therein). We consider here the decay t — > H + b in a model 
with two Higgs doublets. The first order QCD correction to this process has been calculated 
by two groups [|^0| who disagree with each other. In this letter we present a calculation 
of this correction in the limit of a very heavy top quark, i.e., neglecting the masses of the 
bottom quark and the Higgs boson. The effect of finite tuh and mj will be addressed in a 
forthcoming paper. We use dimensional regularization to cope with both the ultraviolet and 
infrared divergences, which greatly simplifies the calculation, especially the real radiation 
part. 

We take Hi and H 2 to be the doublets whose vacuum expectation values respectively 
give masses to the down and up type quarks. The physical charged Higgs H + is a linear 

combination of the charged components of Hi and H 2 , so if we neglect all the Yukawa 

(2) 

couplings except that of H 2 to the third generation = h tt , the top only couples to the H 2 
component of H + . The interaction Lagrangian relevant to the decay t — > H + b is then: 

£ = bg> cos pV tb H + i (^p) b + h.c. 

= ~4 V tb cot (3m t H + i (1 - 75) b + h.c. (1) 

where H + = cos (3H 2 — sin/? if + and cot (3 = (Hi) / (H 2 ) is the ratio of vacuum expectation 
values of the two Higgs doublets. 

In the following calculations we take the space-time dimension to be D = 4 — 2e. The 
mass of the decaying quark is taken to be the renormalization mass scale, and we also use 
it as a unit of energy: m t = 1. In the limit of a very heavy top quark the above interaction 
leads to the tree- level decay rate: 

r(0) (*- H+b ) - %-'J-.rlL 2t r em ^ -ik^ m ^ 2 - (2) 

The first order QCD corrections to this formula arise due to 

virtual gluon exchange and radiation of a real gluon. We first deal with the virtual gluon 
correction to the vertex tH + b. In the limit m b = mn = the spinor structure of this vertex 

2 



remains unchanged and the unrenormalized correction amounts to the multiplication of the 
tree level rate by a factor 



A = C^-^ + - . (3) 



2 r(l + e) / 1 2 
(4tt)t v e 2 e 

where the colour factor CV is 4/3 for £77(3). The counterterm for this vertex involves the 
wave function and mass renormalization constants P.p|: 



If we use the same e to regularize both UV and IR divergencies we obtain for the renor- 
malization constants: 

(4vr)T V e / 

Z>-1 = ^ = 0. (5) 

m b 

The contribution of the virtual correction to the decay rate is: 
T^t (t - H + b) = 2 (A + A c .,) r<°> (t - 77+6 



= 7 =cot a, 3r(2 _ 2e) ^-----12- y j. (6) 

We now turn our attention to the effect of real gluon radiation from the initial or final 
quark. If we denote the amplitudes for these processes by 

Ai and A2 respectively, the contribution of the real radiation to the decay width is: 

r£L (t - H + bG) = ^ cot 2 (3 \V tb \ 2 4na s J dR 3 (t; b, H, G) \A% + A 2 \ 2 (7) 

where the coupling constants have been factored out and t, b, H and G denote the four- 
momenta of the initial and final quarks, charged Higgs boson and the gluon. 

The advantage of using dimensional regularization for the infrared and colinear diver- 
gences is that we need not introduce a mass for the gluon and the integration over three body 
massless phase space is very simple. We choose to parametrise it by the variables x = 2t ■ G 
and z = 1 — 2t ■ b in terms of which the three body phase space integration becomes: 



r 2 4t ~ 7 7r 2t ~ 3 r 1 dr r x dz 

dR 3 (t;b,H,G)= , J , ,J re - (8) 

J K J r (2 - 2e) io (1 - xf Jo z" {x - zf K ' 

After summing over the polarizations of the b quark and the gluon, and averaging over the 
polarizations of the t quark, the squares of the amplitudes become: 

|.Ai| 2 = ~ [2(1 -x)+ x(l -e)(z- x)] 

\M 2 = —(i-O 

x — z 

g 

A l A* 2 + A 2 A\ = — r[l-x + x(l-e)(z-x)] (9) 

2 1 x(x - z) L V /V /J V / 

The integration over the phase space can be done exactly in any dimension. In the limit 
e — > the contribution of the real radiation becomes: 

71 COt ^ 3F (2 - 2e) U 
Although the respective phase space integrations of 



\Ai\ 2 , |^4i| 2 and 

^4i^2 + A2AI give different results from the analogous amplitudes with a W + replacing 
the charged scalar, their sum nevertheless gives the same total contribution in both processes. 
This is in agreement with the general argument based on the equivalence theorem in ref . [|3] . 

Finally we add the effects of the virtual and real gluons to obtain the first order QCD 
correction, so the decay rate (with m t reinstated) becomes: 



Gi 



r (t - H+b) = — ^m? cot 2 (3 \V th \ 



3tt \ 3 



(11) 



This is identical to the result obtained in ref. and also to the analogous correction 
to the decay t — ► W + b 0-§|. If we take a s = 0.1 the first order correction in the limit 
m b — m H — is approximately equal —0.87%. This is in disagreement with the value 
reported in ref. 0,0]. 
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